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the cavity and reservoir energies are equal. A comparison with some particular 
exact solvable cases as, e.g., the Quantum Rem model, is made. 
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1. Introduction 

Quantum Phase Transitions (QPT) arc at the hearth of several fundamental 
phenomena like, e.g.. metal-insulator transition, superconductivity, superfluidity, etc.. 
pp. QPT's can occur both in ordered and disordered models. 

In this paper we show that there exists a very general mechanism for which, in the 
thermodynamic limit, and at zero temperature, the system may undergo a first-order 
QPT through a infinite dilution in the Fock space. Such a phenomena was conjectured 
by C. Presilla and the Author inspired by what they found in a class of exact solvable 
random models [2]. However, these random models are still very particular and do 
not allow to understand what could be the general mechanism leading to first-order 
QPT's piloted by an infinite dilution, if any. 

Let us consider a lattice system of N particles governed by an Hamiltonian 
operator H and let F = {n} be the set of its states, i.e., the Fock space, and H 
the associated matrix representation. The matrix H can always be splitted into two 
parts H = K + V, where V is a diagonal matrix while K connects only different 
states. The matrices V and K play the role of a potential and a kinetic operator 
in the Fock representation, respectively. Let us indicate with \E) and E the ground 
state (GS) and the ground state energy of the system, respectively. Given an arbitrary 
configuration (or state) n, having potential V, let us consider now the reduced Fock 
space IF = F \ n, and let us indicate with \E) and E, the GS and GS energy of this 
reduced system, respectively. Since F and F differ for just one single configuration 
and the Fock dimension |F| grows rapidly with N, we expect that, as the system size 
increases, E and E, will tend to differ for very small terms O(Np), where p = 1/|F|, 
and \E) ~ \E). In this paper we shall prove that this naive result is correct provided 
the potential energy V, associated to the chosen state n, be bigger than the minimum 
possible potential V\. In fact, if V = Vi, when E is smaller than V\, then - as one 
would expect - \E - E\/\E\ = 0(p) and || \E) - \E) || 2 = 0(p), but when E is bigger 
than Vi, once approached the configuration n during a sufficiently long trajectory in 
the Fock space, the system will find energetically more convenient to remain in such 
a configuration for an infinite time. We can see this by using an exact probabilistic 
representation (EPR) of the lattice quantum dynamics [3]. Suppose for simplicity 
that the kinetic term K has the following form K n n i = —r\\ n ^ n i, with 77 > and 
An,n' = 0,1. Suppose also that V be non degenerate. If we define the number of 
active links of the configurations n as A(n) = J2 n ' ^n,n' then, the EPR says that the 
matrix elements of the evolution operator exp(— Ht) 1 t being real, can be evaluated 
as expectations by 



and where n s is the configuration reached at the time s along trajectories extracted 
according to the Markov chain P n ^ n i oc \ n .n'- The steps in the trajectories extracted 
in the Fock space by the Markov chain take place at jump times Poissonian distributed 
with jumping rate rj. This implies in particular that, once reached the configuration 
n, the probability that the system remains for a time t in n, is exp(— riA(n)t). Given 
n, from Eqs. (JTJ) and ^ we see therefore that in general, when t — > 00, up to terms 




(1) 



where the stochastic functional .M' '*' is defined as 



(2) 
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exponentially smaller in t, we can split the expectation and read the GS energy E as 
follows 

Ce~ Et ~Ce- St (l-p) + Ce- Vt p, (3) 

where C, C and C are some finite constants. Now, according to it is always 
E < V\, so that, for t — > oo, the second term in the above equation can be always 
neglected when n is such that V > V\. If instead V = \ \. we see from Eq. (|3j that 
the condition E/N ~ Vi/N for iV — > oo, possibly obtained by varying 77 and the other 
Hamiltonian parameters entering in the matrix V , establishes a critical point through 
which, in the thermodynamic limit, the system undergoes a first-order QPT between 
a normal phase, in which E/N — ► E/N and — ► \E), and a frozen phase where 
instead E/N -> V 1 /N and \E) -y |ra). 

In this paper we will see that, under the condition that p — » faster than 1/N, the 
above simple argument to derive Eq. ((3]) is essentially correct and we can understand 
first-order QPT's as a result of an infinite dilution in the Fock space (not to be confused 
with a dilution of matter) . We will show that such a mechanism is actually much more 
general: given a an infinitesimal portion of the Fock space, the cavity space, and its 
complement, the reservoir space, if - in the thermodynamic limit - the energies of the 
cavity and reservoir space are equal, the system undergoes a first-order QPT. 

We think that this mechanism, due to its wide universality, could explain several 
already known QPT's as well as shed new light toward many interesting less known 
phenomena. 

As an example in this paper we will show how the claimings of our theorem 
coincide with the case of the above mentioned random models [2] , an exact solvable 
class of disordered models which includes, as a particular case, the Quantum Rem 
Model [S]-[7], i.e., a model of disordered spins immersed in a transverse field. 

The main aim of this paper is however to prove the result also in its most general 
form, postponing other possible tests and applications elsewhere. After stating the 
theorem, we will prove it at two different levels: a simple approximate derivation (more 
suitable for the reader unfamiliar with stochastic processes) for the non degenerate 
case, in which finite size corrections are not calculated, and then the exact derivation 
valid for the most general case together with the finite size corrections. 

2. Theorem: QPT's as a result of an infinite dilution in the Fock space 

We now formalize the above phase transition scenario more explicitly. Let F = {n} be 
the Fock space of some system whose physical size (typically the number of particles) 

def 

we parametrize with N and let H be some matrix. In physical systems M = |F| in 
general grows with N as a factorial, however, for what follows, we do not need to 
specify any particular law for M. Let us split H into two parts H = K + V, where V 
is a diagonal matrix while K connects only different states. Let us assume furthermore 
that the Markov chain associated to K (see next section) be ergodic. The matrices V 
and K play the role of a potential and a kinetic operator in the Fock representation, 
respectively |15j . Finally, let us assume that, as occurs in physical systems, both the 

def def 

kinetic terms K(n) = ^2 n , K n n > and the potential levels V(n) = V n . n be extensive 
in the system size N (at least for N — > 00). 

For pedagogical reasons, we find it convenient to separate the theorem in the 
following two parts. 
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2.1. Cavity space associated to the potential level 

Let us order the different levels of V as (Vi < V2 < ■ ■ ■ < V m ), m being the number 
of different levels of V, and let (pi,P2> ■ ■ ■ >Pm) be the associated density of the Vs 
normalized to 1: Yli=iPi = 1- other words 

n 

The number of different levels of V , m, may or may not grow with N. 

Given a possible value of the potential Vi, with 1 < I < m, let us consider the 
following reduced Fock spaces: 

def 

IFi = Un:V(n) = V, n ; ( 5 ) 

which we shall call the cavity space, and its complement 

F= F\U n: v(n)=Vi n ! ( 6 ) 
which we shall call the reservoir space. The matrix H restricted to the spaces F/ and 
F will have some GS and GS energy that we indicate with \Ei) and Ei, and \E) and 
E, respectively. Let us define as e, e, e; and vi, the energy density of E, E, Ei, and 
Vi, in the thermodynamic limit, respectively (for what follows we will always assume 
the existence of these and other similar thermodynamic limits): 

e = lim — , e = lim — , ei = lim —4, vi = lim — \. (7) 

2.1.1. Finite degenerate case Let us suppose first that the level Vi be finitely 
degenerate in the limit N — » 00. Note that this implies in particular that be — > 
for N — > 00. Let be p ; 7V — > 0, then for / > 1 

e = e, (8) 

while for / = 1 

' ^' ^L- Wl ' (9) 
ui, e > v\. 

Eqs. ([H]) and © tell us that, in the thermodynamic limit, the GS of the reduced 
system approaches the GS of the full system if I > 1, while if / = 1, we have that the 
equation 

e = «i, (10) 

as an equation in the space of the "Hamiltonian" parameters entering H , represents a 
critical condition where a first-order QPT takes place: we will prove that, for e < v\ 
- as in the normal case I > 1 - the GS's \E) and \E) differ for terms O(pi) 1 ^ 2 , and 
\E — E\ — O(Npi), while for e > ui the system gets frozen and \E) tends to n±, the 
configuration associated to the lowest potential level V\ . 



2.1.2. Infinite degenerate case Let us consider now a more general case in which the 
level Vi may be infinitely degenerate in the limit N — > 00. Now, the cavity space 
associated to the level Vi has, not only an extensive potential energy, but also an 
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extensive "kinetic energy". Let the states in F; be connected. Then, if for N — > oo, 
piN — > 0, we have 



For 1 = 1, Eq. pip is analogous to Eq. ^ and represents the natural generalization 
of the former case to the infinitely degenerate case, the critical condition now being 



Note however, that unlike the former case, we cannot assume the analogous of Eq. 
l[5|) for I > 1, and in particular here we are not able to know a priori whether a more 
general critical condition e = ei may or may not have solutions. However, whatever 
be I, according to Eq. (fTTjl . if the equation e = e; admits solutions, the system will 
undergo a first-order QPT phase transition with the same mechanism as in the finite 
degenerate case with a normal phase and a frozen one where the GS lives in the Fock 
space Fj. 

2.2. General case 

Finally, we have the following natural generalization. Let F be any connected subspace 
of F having a fraction p of configurations with pN — ► for N — > oo [TB] , then in the 
thermodynamic limit, for the energy density e we have 



where the density energies e and e are those associated to the spaces F and F, 
respectively. As in the previous case, we are not able to say nothing about the 
possibility that the equation e = e has or has not solutions in the space of the 
Hamiltonian parameters present in the original matrix H. Again, in the affirmative 
case, we have a first-order QPT taking place with the same mechanism. 

Concerning the finite system in which then p is small but finite, we will also 
prove that, near the critical point, the difference between E and E, or E, is of the 
order 0(^K out ), K ont being the kinetic energy by which, once entered in the cavity, 
with asymptotic probability 7f, the system is expelled out from the cavity space itself 
(i.e., the space F). More precisely, ff is the invariant measure of the Markov chain 
associated to the matrix elements of the "kinetic" matrix K. Note that, as we will 
point out later, since under the hypothesis of the theorem the matrix K is connected 
and the kinetic terms extensive, it is always ff < p, for N — > oo, so that the condition 
pN — > ensures also that be nN — > [TS]. The detailed expression for K out is quite 
cumbersome and is left in the proof (see Eq. ([71]) ). The energy K out , besides to be 
dependent on the kinetic terms associated to the boundary of the cavity, depends also 
on the overlaps that the reservoir and cavity GS's have on their internal boundaries d¥ 
and 9F, receptively The larger are these Fock-boundaries contributions, the stronger 
will be the finite size effects. 

3. Proof 

3.1. Exact probabilistic representation (EPR) 

The derivation of Eqs. (O, © and (fT5|) follows from an exact probabilistic 
representation derived in [3] and that here we recall briefly. In general, given an 




(11) 



e = 



(12) 




(13) 
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arbitrary Hamiltonian matrix H , each row n of the corresponding hopping matrix K 
has a different number of non zero elements, let say A(n). We call A(n) the number 
of active links of the state n (in graph theory A(n) is called degree of the vertex n). 
Given the matrix H, and its hopping and potential matrices, K (the non diagonal 
part of H) and V (the diagonal part of H), respectively, we define a virtual dynamics 
as follows. Let us parametrize the matrix K as 

such that \\n,n>\ can be either or 1 and rj n _ n ' > 0. In graph theory the matrix with 
elements |A n r »'| is known as the adjacency matrix: in fact, it establishes whether two 
given states n, n' are first neighbors or not. We consider the Markov chain defined by 
the transition matrix P with elements 

p„,„. - !^!, (i5) 

where 

R(n) =]T \K n , n , | . (16) 

to' 

Starting from a given initial configuration no, we draw a new configuration n\ with 
probability P nQ , ni - By iterating this procedure for TV steps we construct a path, 
or trajectory, in the space F no,rii, . . . , rajy- Note that, since P n ,n' is supposed to 
be ergodic, after an infinite number of jumps the Markov chain converges to the 
invariant measure 7r, which coincides with the left eigenvector of the transition matrix, 
7r T P = 7r T . It is simple to verify that [4] 

Notice that, the hypothesis that the Markov chain be ergodic, implies in particular 
that for any n be R{n) > 0, so that from Eq. (|TT|) we have 

, , R(n) 1 , . 

7T (n) < ^ , - T-r, 18) 

V j - mm n , eF R(n>) |F|' V ' 

and more in general, if F is a subspace with a measure p coming from Eq. (0]), we 
have 

- rfe, f = Sn e r^(") < max„ sF fl(n) 

En R ( n ) ~ niin„ eF i?(n) * V ; 

so that, if the i?'s arc all extensive, and if, in the thermodynamic limit, p — > 0, also 
7f — > (see note [H]). 

Along each finite path with N steps we have the sequences of data Aq, A\, . . . , An, 
Ro,R!,...,R N , V ,Vi,...,V N , Ai,...,Ajv, and r)i,...,r] N , where 

A k =A(n k ), k = 0,...,N, (20) 

R k = R(n k ), k = 0,...,N, (21) 

V k = V(n k ), k = 0,...,N, (22) 

Afe = A TlfcliTlfc , k=l,...,N, (23) 

»/fc = »? n *_ ll n Jb , k = l,...,N. (24) 

Note that, by definition, along a path of the Markov chain |Afc| = 1. 

If at each step of the Markov chain we extract randomly a jump time s 
distributed according to a Poisson process of parameter p, and along each trajectory 
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we indicate the sequence of these drawn jump times as si, S2, . . . , Sjv t) Nt being the 
number of jumps occurred along the trajectory, then the following exact probabilistic 
representation holds [3] 

(n\e- m \n ) = E(8 n , nNt M^), (25) 
where the stochastic functional AdnlP is defined as 

M^'V = \ T[ —Xke [pAk - 1 ~ Vk - 1 ^ Sk ~ Sk - l)> \ e l pAN *~ VN *'l( t ~* N *' ) (26) 

if N t > 0, and Mnf = e^ pA °~ Vo ^ if N t = 0, and the symbol E(-) means expectation 
over the Poisson processes and the Markov chain. In Eq. (|26|) . we put so = 0. 
According to this representation, we have in particular 

Y / (n\e- Ht \n ) = E(Mt t) ), (27) 

n 

so that we can evaluate the ground-state energy as 

E = lim -9 t logE(^' °)). (28) 

t—*oo \ u / 

Iii the following we rewrite Eq. (|26p in a more compact form as 



(29) 



More in general, if at each active link (n, n') we associate a Poisson process with rate 
pn.n' dependent on the link, the following generalization holds: 

m£v = ( e m^s)-v{ ns )]ds tt m x \ (30) 
V fe=V fc / 

where R(n), the weighted number of active links, is defined in Eq. (|16|) . For what 
follows we will always make the choice 

Pn.n' = Vn.n' (31) 

so that the product in Eq. IpTl]) simplifies in 

M [o,t) = e ;«[JK»,)-v(».)]i. -q Xk (32) 

fc=i 



3.2. Proof of Eqs. {Bp-flfy (simple derivation) 

Let us consider the degenerate case, i.e., the case in which the cavity consists of a 
single state n. In Sec. II we have indicated with E the GS energy of the system 
living in the reduced Fock space F = F \ n and having (for n ^ n) the same matrix 
Hamiltonian representation H (alternatively, such a system could be still seen as a 
system living in the original Fock space F having as operators K = K and V n n = V n ,n 
for n ^ n and V ni n = +00). This Hamiltonian in the following will be indicated by 
H. We find it convenient to formally define also an Hamiltonian associated to the 
single state n and taking (for the moment being) the single value H = V. 
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Given an initial configuration no G F [TH], let fc< be the stochastic variable 
counting how many times the system, starting from no, enters in the state n during 
the time interval [0,t). We can decompose the expectation as 

oc 

E (m^) = E (M^\k t = k) Q t (k), (33) 

fc=0 

where in the rhs we have introduced Q t (k), the probability that the system enters 
in the state n exactly k times during the time interval [0, t) (dependencies on the 
initial configuration are understood). Let t\, t%, . . . , tk be the jump times at which 
the system enters in the state n through the jump of some active link having signs 

An _ .n + , ■ ■ ■ , A„ _ , n , , respectively, and t±, T2, . ■ . , the corresponding living times 

H H _ *fc ' *fc 

in n, from which the system leaves through a further jump of some active link having 
sign \ n n A n n _i_ , respectively. Here, to avoid confusion, we 

have introduced the symbols for the configurations immediately before the jump taking 
place at the time tk, as n t - , and immediately after the jump as n t + , and similarly for 

k k 

n t , T - and n t +T +. Then, from Eq. ((321) we see that Eq. ([53^1 becomes 



k=0 



f l+ T 



xM 



tfc-l+Tfc— l,tfe 



A Tl _ I „ + ^ [t *'* k+7 * ) A n _,„ + X^+^^|fc t = fc)g f (fc), (34) 

where the stochastic functionals .Ml '*- 1 and A^' '*-* are defined similarly to the 
functional yVf[°'*) introduced at the beginning but the Hamiltonian H is replaced with 
the Hamiltonian's H and H, respectively. Let us note that the two set of stochastic 
variables ti,t2, ■ ■ ■ , tk and ti, T2, . . . , Tk are independent. Furthermore, the variables 
T\,T2, ■ ■ ■ ,Tk are almost each other independent, the only constraint being that their 
sum cannot overcome t. So that, taking into account that once entered in n, the 
infinitesimal probability that the system leaves n after a time between r and r + dr 

is dTRexp(-RT), with R d = i?(n) [17], according to Eq. (1321) we have 

°° /t 

E 



(Xl -*)) = V / d Xl dx 2 ...dx k E (m^^M^uU) . . . M[tk+* k ,t) 



where J dx\dx2 ■ ■ ■ dxk is a shorthand notation for the integral extended to the only 
fc-dimcnsional set where the variables satisfy the constraint < x\ + x 2 + . . . + Xk < t, 
and we have introduced the following stochastic normalized average 

k \ 

def 




\i=l 

E fM[o,ti)Mlti+xi,t 2 ) . . . M\tk+x k ,t)\k t = kj . 

Note that, in the above expression, we have for brevity written the expectations 
appearing in the numerator and denominator as non conditioned by the times 



(36) 
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Ti = xi, . . . ,Tk — Xk- However, in the limit t — > oo, the stochastic normalized average, 
which is the ratio of the expectations, becomes a constant independent on the t's. 

We observe that, in each trajectory, the total time t during which the system is 
not in the state n is t — t — (xi + . . . + Xk)- Furthermore, for any fixed value of k, each 
interval, [ij + xi, i;+i), in which the system is not in the state n, has a length which is 
of the order 0{t/k). Therefore, when t — ► oo, the stochastic functionals A^' t,+:E! '* ,+1 - 1 
in Eq. (f3l>|) get uncorrelated and for any finite k we are left with [T§] 

£^M [aM) M [tl+Xut2) ---M tk+Xk ^\k t = k;n =x u ...,Tk = x^j 

^ C t (k)e-^ t - (xi+ - +x ^\ (37) 

where E is the GS energy of the system with Hamiltonian H, and Ct(k) is a bounded 
quantity smoothly dependent on t. In the proof for the general case we will show that 
both these constants and the stochastic averages (f3l>| are well bounded functions. 
More precisely we will show that the net result is given by 

C t (k)R k (j[\ tl \ tl+x }j =C no {-K ont f{l + o{e-^-^ t )), (38) 

where: E\ > E is the energy of the first excited state of the reservoir, C no = O(l) is 
a constant depending only on the initial condition no, and K ont , has the meaning 
of a kinetic energy pushing the system out of the configuration n. It is in fact 
K \a ~ — (R + — R~), where R + and R~ are the kinetic terms obtained as sums 
over the positive (A = 1) and negative (A = — 1) active links, respectively. 

By using now Eqs. (|37|) and (j38j) in Eq. (|35|) . up to terms exponentially smaller 
in t, we get 

/ \ ~ - 00 r' 1 ~ - 

E{M^) =C no e- Et J2(-K out ) k dx 1 dx 2 ...dxke^ E - v ^ + - +x ^Q t (k). (39) 

We have now to understand the behavior of Qt(k). We recall that the total 
number of jumps made by the system during the interval [0,t) is of the order (R)t, 
where (R) is the average of the weighted number of active links R along a trajectory 
making an infinite number of jumps [20] . Now, the way the system jumps from one 
configuration n to the next one n' is a finite ergodic Markov chain Pn, n '- Therefore, 
there exists a critical value N c for which, if the number of jumps N satisfies N 3> N c , 
up to terms exponentially smaller in N one has (P N ) n ,n' = t(w'), where 7r is the 
invariant measure given by Eq. (|17p . In particular, this implies that, for N 3> kN c , 
up to terms exponentially smaller in N, the probability VN{k) that the Markov chain 
visits k times the state n is pAr(fc) = (7r(n)) fc . Therefore, there exists a critical value 
t c ~ N C /(R) for which, if t 3> kt c , up to terms exponentially smaller in t, we have 
also Qt(k) = (7r(fi)) fc [3T|. Furthermore, recalling that < Qt(k) < 1, and that 

/ * dxi . . . dxk = t k /kl, we see immediately that the series in the rhs of Eq. is 
absolutely convergent for any t. This fact allows us to pass the limit of t large under 
the series symbol and to analyze the behavior of the system by simply analyzing term 
by term the single contributions of the series. The integrals involved in the rhs of 
Eq. (|39p cannot be expressed in a closed form. However, the analysis of the first few 
already makes clear the general structures. For V ^ E we have: 

E (M^) = e- Et C no Q t (0) + (e- Vt - e~ St ) j^C no (-K out ) Q t (l) 
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e~ vt t- 



E-V 



,-vt 



-El 



E-V 



C no (-Xout) Qt(2) 



+ • • . + C no £ fe (t) (-JT out f Q t (fc) + . . . , (40) 

where Qk (t) is a multivariate polynomial in the variables exp(— Vt), exp(— Et), and 
t,t 2 . . . , t 1 . Let k t be an integer such that 1 < k t <C i/i c . By using Qt(k) ~ 7r(n) fc , 
for 1 < k < kt, and then Qt(0) ~ 1 — 7r ("') fe f° r t large enough, up to terms 

exponentially smaller in t, Eq. (|40|) becomes 

E (.M^)) = e- St \~_ 2 ^ C no + (e~ Vt - e^) ^-rC no (-K out ) ir(n) 



E-V 



e~ v H- 



1 



■,-vt 



-,-Et 



E-V 



1 



V 



+ ■■■ + C no G k (t) (-K out ir{n)) k + ... + C no G- kt (t) (-ET out ^(n)) £f + K u 



(41) 



where IZt is the rest of the series taking into account the contributions for k > kt- 
Note that, by construction, for t — > oo, k t — ► oo, 72^ ^0 and Eq. (|4T|) becomes exact. 
Concerning the case V = E instead we have 

E (M^A - -~ Vt ( 



I — 2?r / n ? C f no + iC™,, (-if ut) t(ti) + %rC no (- K out n{n) f 
1 — 7r(n) Z 



Cno Y& 



t 7r(n)) 



+ . 



(-.Kout7r(n)) 



(42) 



and must coincide with the analytic continuation of Eq. (|4Tj) for y E 1 . According 
to Eq. (|28|) . we have now to analyze the behavior of the rhs of Eqs. (|4Tj) and |42|) in 
the limit t — > oo (note that, if we want to calculate the GS, we are allowed to take the 
thermodynamic limit only after performing the limit t — > oo). Were the series involved 
in these equations truncated at a finite number of fc, we could get immediately \E) 
and E, by simply choosing between \E) and -E, or |n) and y, according to which is 
the smallest one between E ov V . This would erroneously bring us to claim that the 
phase transition scenario summarized in Sec. II occurs actually also for finite sizes! 
Such an error of course depends on the fact that we have truncated the scries. In fact, 
when we have to sum an infinite number of terms of the form exp(akt), in general, 
the resulting sum will be a new exponential in t but having a new suitable exponent. 
We can better understand this phenomena from the easier case V = E, i.e., near 
the critical point. In this case we are able to sum the series in Eq. (|4"2"|) and we get 
E = V + TrK out and not E = V . We see therefore that: i) only in the limit tt — > 
we are allowed to claim from Eqs. |4T|) and (|42| which is the GS and the GS energy 
by simply taking the smallest between E or V, obtaining then the phase transition 
scenario described in Sec. II; ii) there is no phase separation for 7f finite; Hi) at the 
critical point, or in a sufficiently small region near the critical point, for small but 
finite n, the GS is always a mixed phase to which contribute both the reservoir and 
the cavity, and E — E ~ E — V is ot the order 0(^K out ), so that, given n, the bigger 
is -fTout, the stronger will be the finite size effects. Note finally that, due to Eq. (fl"8]) . 
all these statements can be reformulated in terms of of p = 1/|F|, being pN — > a 
sufficient condition for having nN — > 0. 
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3.3. Proof of Eqs. fTI \) -l[T3 )) (exact derivation) 

We now generalize the above derivation to the case in which the cavity F is any 
infinitesimal subset of the Fock space such that ttN — > for N — > oo. It is convenient 
to define the internal boundaries of the reservoir and cavity sets as 

<9F = U ne f: 3n ^f : K n>n ,±o n - ( 43 ) 

^ = J U n£?:3n , #: ^ x#0 n. (44) 

For each configuration n £ F we decompose the total number of active links A(n) as 

A(n) = A in (n) + A out (n), (45) 

where Aj n (n) and A out (n) represent the number of connections of n to configurations 
belonging or not belonging to F, respectively. More explicitly 

AM d = E (46) 

n'SF 

Aout(n) = f lVn'1- (47) 

Clearly, if n 6 F \ <9F then A out (n) = 0. Similarly, for n 6 F, we decompose the 
weighted number (or kinetic term) R(n) as 

i?(n) =i? in (n)+i? out (n), (48) 



where 



i? in (n) £ (49) 

Rout(n)= f J2 \ K n,n'l (50) 



The corresponding definitions for the reservoir space are understood. Sometimes, to 
stress whether the kinetic terms i?; n and i? ut refer to the cavity, or to the reservoir, 
we will write i?; n and i? ut, or R\ n and i? ut, respectively. Note that the correct kinetic 
terms associated to the Hamiltonian's H and H of the cavity and reservoir spaces, are 
-Ri n (n) and i?i n (n), respectively. 

Let ti , . . . , tfc be the stochastic times at which the system enters in the subspacc 
F, and let n,. ..,Tfc be the corresponding living times. Note that, by construction, 
at the times tj", . . . , ijj", as well as at the times ti + , . . . , t^ + r^T , the system is 
on the boundary d¥. As in the previous paragraph we can repeat the same steps 
also for the present degenerate case, Eqs. (pt5| and being formally identical with 
the obvious generalizations for the symbols Qt(k), A^I '*-* and M\ 0,t \ Note however 
that .A/f' '*' and M\°^\ defined as restrictions of the stochastic functional M.^ '^ to 
the configurations of the spaces F and F, respectively, are not yet associated to the 
Hamiltonian's H and H defined as the restriction of H to these spaces. In fact, 
the kinetic terms i?'s appearing in these stochastic functionals arc different from the 
kinetic terms i?i n 's associated to H or H when the system is on the boundary of F 
(<9F) or on the boundary of F (<9F), respectively. Furthermore here, to repeat the 
analogous derivation of the previous paragraph, we should calculate exactly the living 
times in the subspace F with respect to which we have explicitcd the expectation (i.e., 
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the analogous of Eq. |55|) ). Let us first consider the simplest case in which for any 
n G F we have a fixed value for both the potential, V(n) = V, and for R{n) = R, 
i? in (n) = -R; n , and then also i? ut(«) = -Rout- We observe now that, once at the time 
r = the system enters in a configuration n G F, the infinitesimal probability that 
the system leaves F after a time between r and r + dr is given by drR ou t exp(— R out T). 
Therefore, according to Eqs. ([52"]) and we see that the analogous of Eq. (|35p is 
given by (for simplicity in this example we consider a bosonic case so that the A's 
associated to the jumping links are all positive) 

E[M l °^) =Y\I d Xl dx 2 ...dx k E(M lQM) M ltl+Xl - t2) ■ ■ ■ M [tk+Xk > t] 

x e [-^+v]t*^-+^)\ kt = k;n = xi, . . . ,n = x k ) 

x R k out Qt(k). (51) 

We see then that, similarly to the non degenerate case, the probabilities of the living 
times in F cancel exactly with the portions of the kinetic terms -Rout, while the rest 
of the kinetic terms, i.e. i? nl , are exactly those associated to the true H (note that in 
this example E = — R- m + V). Now, for the general case, analogously to the previous 
paragraph, we can write (we suppose to start from a configuration in F) 

OO / 

x M^ k+Tk) X + + + Ml k+Tk f\k t = k) Q t (k), (52) 

We would like now to generalize Eq. (|5ip . However, the probability of the living times 
in F for the general case represents a formidable task. Nevertheless, as we shall see, 
for our aim we do not need to calculate this quantity, but only a suitable conditional 
expectation. In fact we can use the following Lemma. 
Lemma. 

Let t be a generic random time at which the system, for the first time before t, leaves 
some region of the Fock space (for example r can be the first time at which the system 
leaves, or enters, the cavity). Then we have 

E(M& T) ) = f o dxE(^R out {n x -)\r = x) , (53) 

where the symbol M^.^ indicates that in the stochastic functional now the kinetic 
terms are only those internal to the region: 

M[°£ o d M ^[RM-vi^ds -Q Xk _ (54) 

fe=i 

More in general, for what we need, we have 

\ n ° H ■ *t t+ tl+T! h+t+ 

X M [ tt" +Tfe) An _ ,„ + M [ l k+Tk : t] \k t =k] = 

' dx\ . . . dxk E (^££>* tr n tt M [ £^ ■ ■ ■ 



Quantum Phase Transitions induced by Infinite Dilution in the Fock Space 13 
x \kt = k;r 1 =x 1 ,..., r k = x k ) . (55) 



Furthermore we have the following crgodic property 



lim V _ ' ° rr '- = 1, (56) 



y l m;n J 

and similarly 



lim V ,~ rn .,^ ' =1- (57) 



Vi— »o° 



\ m;n 



The proof of this Lemma is reported in Appendix. 

Let us now exploit the Lemma for the case k t = 1 [22]. By using the Markov 
property we can rewrite this contribution as [5] 

x^t£« Bl '* ) + l*t = 1 ;n = a J i > )= (58) 
xE fe^Jn^ |fc = l;n = ^ . (59) 



m:n 4. 1 

The last conditional expectation that appears in the rhs of the above expression is a 
random variable that, using the EPR, can be written as 

E (^r 1 ^ K-h* ) = ^ ti+ ^e-^-^+-)l + O (e-*I«-^M) , (60) 
where E\ is the first excited state of the reservoir, and is given by 

C n = f ^2(ri\E)(E\n). (61) 

n' 

Note that, since \E) is normalized to 1, it is Ca = 0(1). By inserting Eq. (pOj) in 
Eq. (|58p. up to terms exponentially smaller in t — (ti + xi), and by using again the 
Markov property, we get 

E (^K;^ l) ^t(n tl+Br )A nti+B _, nti+Ti+ C Rti+- Jn t +) |fc t = l;n = a*) . (62) 
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In turn, the last expression can be expressed in terms of the probability that given 
n t 1 +x~ > tne system goes from n ti+x - to a generic configuration n 6 F. Taking into 
account Eqs. (TT4")) - (rTrj)) , we have 

E r u * ? — — ^^j- 

E (vW|*^+ Xl Vt+;n = *l) (E " Jif « tl+ . r .» ^cavity, (63) 

where, in the last expression, we have introduced the normalized average 

^ + - ■ ri C'n) cavity — 



(64) 



Note, in Eq. (|63|) . the presence of the denominator i? ut- This is due to the fact that, 
the probability to be used in Eq. (|63[) must take into account that the expectation 
to be calculated is conditioned by t\ = X\, i.e., the system at the time x\ leaves 
the cavity with probability 1, and, with a jump, can go to any of the R ut(n tl+X -) 
configurations of the reservoir space (in the non degenerate case -Rout = R)- 

Note that, as x\ goes to oo, the normalized average (|64[) stabilizes exponentially 
fast to a asymptotic value independent on the initial condition n t +. Therefore, taking 
into account that eventually we need only to keep track of the behavior for t — > oo, 
and then also x\ —> oo, up to exponentially smaller terms, the rhs of Eq. (|62p becomes 
now 

E (e- £[t - (il+ - )] A?!S;An tr ,n 4 E (^+ Sl) K+;n = *x) \h = l;n = Kl ) 

lim (} -K n n C n ) cavity (65) 



X 

t — >oo 



Similarly to what we have done before in Eq. (|60|) . we see that the last conditional 
expectation that appears in the above expression is a random variable that, using 
again the EPR and Eq. (J^H) , can be written as 

E (^fc^'l^n = xx) = C^e-^ + O (e~^) , (66) 

where E\ is the first excited state of the cavity, and C n is given by 

C n = f Y^( n '\E)(E\n). (67) 

n' 

Note that, since \E) is normalized to 1, it is C n = O (1)- By inserting Eq. (fl)!)]) in Eq. 
(1651) we obtain 



E (e- S[t -^ ] Mf n ^X nt _^ +E (^fe^K+;n = xx ) I*. = IjTi = a* 
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Finally, by performing analogous steps as those done in Eqs. (|63|) -(|65 p . and by using 
Eq. ([57]) . we arrive at 



E(e-^-^ + -)l^-AlSA„_,„ + C ftt+ |fc t = l;r 1 =x 1 ) = C no x 



1 1 

-a; 



E f e -Elt-(t 1+ ^-E Xl -Et llkt = 1 \ Hm { J- _ C„) rcBorTOir , (69) 

where we have introduced the normalized average (see note |22j ) 

v E ( Jw! '?. V - "" r '" C*„|n ; fc t = l:t x = j/x | 
-Jf n _,n , , I m;rl ° ^ n flout (n ) 711 u ' 11 ' 7i I 
/ V "i r \ . __\ !i /_ fyn^i 

\7 J ~ . . ^n/ reservoir / _ m \ \ 'V / 

n J2o«t(n Vi - ) E |n ; fc t = 1; = Vl J 

Note, in Eq. f?0]) . the presence of the denominator i? ut- This is due to the fact that, 
the probability to be used in Eq. (|69|) . must take into account that the expectation 
to be calculated is conditioned by t\ = yi, i.e., the system at the time y\ leaves 
the reservoir with probability 1, and, with a jump, can go to any of the R ut{n t -) 

configurations of the cavity space (in the non degenerate case i? ou t(^ t -) = In- 
putting all together from Eq. ([55]) to Eq. ^ we hence obtain 



E^!::^„ rn ^!^\ ut (n (l+ir )A„ i+irIlti+r , 

xMtt X1 ' t \\h = l;r 1 = x 1 ) =C no e~^-^ (71) 



m;n + ■ 



E~ ^ — / V — 
C n ) cavity Hm ( / j — ~C n ) 
reservoir- 
(72) 

n n n out (n t -) 

Notice, at this leading order, the null role played by the random variable t±. 

In general, by iterating A; times the above procedure, up to terms exponentially 
smaller in the random variables t and {x\ + . . . + Xk), we get 

[0,*i) \ d J„ \\ j^[tk,t k +x k ) 



nv.n 



ft 



x Rout(n t + -)X n _,« ++ Mf n k ~ Xk - t] [\k t =k;n =x 1 ,...,n = x k 

C no e- St+ ( S - S ^ + - +X "\-K out ) k , (73) 
where 

d& ^ ^ — ' ' — ^/\. j *W — 

K OVLt = ~ Hm ( > K n _,»C«)cavity lim (> — '"' - C n ) reservoir- (74) 

t— *oo * — ' * i^oo * — » f? , in \ 

„ „ -rtout l«t-J 

We can now insert Eqs. (J73]) and (J7I|) in Eq. ([55|). Finally, from Eq. (JSH), up to 
terms exponentially smaller in t, we arrive at the analogous of Eq. p9|) given now by 

E (^[°'*n =C no e- Et 2~2 dx 1 dx 2 ...dx k e^- E ^ + - +x "H-K out ) k Q t (k). (75) 



We see that Eq. (|75j) is just the wanted generalization of the specific case we had 
analyzed before in which for any n € F we had a fixed number for i?(n) = R, 
Rin(n) = R in , and R ut(n) = R out (see Eq. (JSTJ). 
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The phase transition scenario described in the theorem by Eqs. ()lip -(fT5 |) follows 
with the same identical argument given in the previous paragraph (the non degenerate 
case). We see in particular that, for a finite system and then a finite value ff, near the 
critical point, the difference between the energy of the full system E, and the energy of 
the reduced systems E or E, is of the order 0(K out 7r). We observe that, using the fact 
that the GS's \E) and \E) are normalized to 1, from Eq. (|74[) . we have approximately 

K(mt ~ — (-Rout — -^out ) cavity j (76) 

where R^ ut and Ra Ut , are the kinetic terms obtained as weighted sums over the positive 
and negative active out-links, respectively Since, in general, —((R + — R~}) is just the 
kinetic energy of the system [5] , we see that the meaning of K out is that of the kinetic 
energy by which the system, once entered in the cavity with an asymptotic probability 
7f, is expelled from the cavity itself. Given a small but finite value 7f, the lower is this 
energy, the greater will be the ability of the system to remain in the cavity for a 
longer time, making then faster (as TV grows) the approach to the phase transition. 
Note also that, according to their definitions, the constants C n and C n , given in Eqs. 
(jSTj) and (jCTTf . measure the overlap of the GS's \E) and \E) with the reservoir and 
cavity internal boundaries df and OF, respectively. Since, roughly speaking, these two 
overlaps appear as a product in -Kout, we see that, as is intuitively evident, the larger 
are these overlaps, the stronger will be the finite size effects and, more precisely, given 
the sum of the two overlaps, their effect will be maximal when both the overlaps are 
relatively large, while minimal when one of them is relatively small. 



4. An exact solvable example: The Random Potential Model 

In [2] we have solved explicitly a model, called the Random Potential Model, a 
quantum quenched-disordered model, in which the potential levels Vi, with I = 
l,...,m, take independent random values distributed according to an arbitrary 
assigned distribution pi, whereas the kinetic matrix K is arbitrary with the only 
constraint that its associated Markov chain be ergodic. We have found that, up to 
small corrections that become negligible in the thermodynamic limit, the GS energy 
E of this model satisfies: 

where E^ is the GS energy of the system with zero potential (V = 0), and the 
constraint E < Vi makes the solution unique. We have then analyzed in detail this 
equation in the special limit p\ — > and we have seen that the system undergoes a QPT 
first-order phase transition in correspondence of a suitable choice of the Hamiltonian 
parameters (i.e., the values of the Vi, their weights pi and the parameters entering the 
matrix elements of K that affect then E^- '). Let us consider the simplest non trivial 
example with only two potential values V± and V2. In this case, Eq. (|77p is a quadratic 
equation for e = E/N which can be solved explicitly. Observing that pi = 1— pi, and 
defining also e^ - 1 = E^ /N, v\ = Vi/N, and V2 — V2/N, and taking into account that 
e^ ) is negative, one has [5] 
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the other solution being incompatible with the condition eo < V\. We see that, for 
any p\ > 0, eo is an analytic function of V2 — v\ and e*- -* . However, for p\ — > we have 



lim eo — f i 



(v 2 -Vi + e (0) ) - |i> 2 - Ui + e (0) 
V2 — Vi + e^°\ V2 — V\ < — e*- ) 



(79) 

I 0, v 2 - vi > -e (0) , 

i.e., in the space of the Hamiltonian parameters V2 — v\ and e' -*, a singularity shows 
up at the critical point 

V2 - Vl = - e (°), (80) 

On the other hand, we see that Eq. ([80]) can be written also as 

e = VL (81) 

More in general in [3J, we have shown that, for p± — ► 0, the critical condition can be 
written as 

w = Wv m 

where the function W is defined as 

m 

W(p 2 ,V 2 - Vi,... ,P m ,V m - V\) = V" , (83) 

and vi = Vi/N. Again, from Eqs. (|77|) and (|83|) . and taking into account that, by 



definition, p\ = in the reservoir space, we see that the critical condition (|82[) is 
nothing else but the equality of the density energies e and v\ that we have defined in 
our theorem for the finitely degenerate case. 

An important application of this model is the Quantum Rem model, i.e., the 
Dcrrida's model [9] immersed in a transverse field, and the QPT here signals the 
transition of the system between a quantum paramagnetic state and a glassy one [5]- 
[7|. These kind of models have important connections to the general problem of the 
Quantum Annealing [TU] . 

For a general model, having then non random potential levels, Eq. (|77|) is of 
course not valid. However, it is possible to develop a perturbative approach whose 
zero level approximation coincides with Eq. (|77|) , while non trivial correlations of the 
potentials and the kinetic terms are taken into account at higher order of the theory 
[11] . Since the general structure of the equations of this perturbative theory remains 
similar to Eq. ((77)) itself, we believe that a phase transition scenario similar to the 
one above described will take place in the infinite dilution limit of suitable regions of 
the Fock space. 



5. A comment on the first exit times and the EPR 



Concerning the EPR and the random times r's, as times at which the system leaves 
the cavity for the first time before t, a comment is in order. From the EPR we 
know that to calculate the GS energy of the system, we need to know the behavior 
of the expectation of the stochastic functional involved, only in the limit t — > oo. 
Nevertheless, in proving our theorem exactly, we needed a perfect balance between 
weights and corresponding probabilities at each finite time interval. For example, from 
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the formulae present in Appendix, it is not difficult to recognize that, when t — > oo, 
the pdf p(r) for the random variable r to leave the cavity for the first time before t, 
satisfies the following asymptotic law 

lim l0g[p(T)] = -ET, (84) 

r — >oc T 

where E* is the GS energy of a "bosonic" Hamiltonian restricted to the cavity and 

- def 

defined as follows: H* = K* + V* , where, for any n,n' G F, if* n , = \K n ^ n i\, K n ^ n > 

def 

being the original kinetic matrix of the given system, and V*(n) = R(n), R(n) being 
the (total) kinetic terms associated to K n n i. Note that, by construction, the kinetic 
terms of H* are given by R*(n) = R\ n (n), whereas R(n) = R m {n) + R out (n), with 
i?i n (n) and i? out (n) defined in Eqs. (|49|) -([50 ]l . Notice that: i) if in the place of 
V*(n) = R(n) we had V**(n) = i? in (n), the GS of the corresponding H** would be 
the uniform state with GS energy E** = 0, ii) for the same reason, and due to the 
fact that the cavity space is connected to the reservoir (one of the hypothesis of our 
theorem), it is E* > 0. Interestingly, properties i) and ii) can be seen as classical 
results of algebraic graph theory [T2], but they can also be derived immediately by 
using the EPR [23] . We see here clearly the connection between the pdf p(r) and the 
condition E* > 0. The fact that the cavity be connected to the reservoir ensures that 
the probability to escape from the cavity in a finite time be finite, and that the decay 
of the pdf be exponential. 

Now, the pdf given in Eq. ([54"]) is an interesting element in statistical mechanics. 
Nevertheless, it is not useful for our theorem: in our theorem we need a less difficult - 
but still exact - pdf as the p(r = x) given in Eq. (|A.5[) . The perfect balance between 
weights and probabilities, condensed in Eq. (|53"|) , is a feature of the EPR that other 
approximate probabilistic representations, based, e.g., in Trotter-Suzuki formula (see 
e.g., [13]; see also [12): cannot satisfy, and a tentative to prove the theorem with 
them would fail. The exactness of the EPR is the crucial ingredient at the base of the 
first-order QPT's scenario stated in the theorem [24] . 

6. Conclusions 

Starting form an exact probabilistic representation of an arbitrary quantum dynamics 
[3J, a very general theorem on QPT phase transitions has been here derived. The 
critical equation, which consists in the equality of the reservoir and cavity energies 
E and E in the thermodynamic limit, under the condition that the cavity space F 
be infinitely diluted with weight p = |F|/|F| going to zero faster than 1/N, has a 
very appealing physical explanation as follows. Once the system enters in the cavity- 
subspace F, if its energy E is lower than the energy of the reservoir-subspace E, just 
due to the fact that p — > 0, the system prefers not to leave the cavity anymore, since, 
otherwise, once left the cavity, the probability for the system to find again the cavity 
would be very low and, during its traveling in the Fock space, it would accumulate 
an higher energy E. On the other hand, when instead E < E, the system is no more 
forced to be in such a rare event in the Fock space, so that, just due to the fact that 
p — > 0, the cavity is almost not visited and the system remains in the normal phase 
with an energy given by E ~ E with a relative error of the order 0{p)- 

We point out that this phase transition scenario picture takes place in the Fock 
space and that the required infinite dilution pN — > has nothing to share with a 
material dilution. In fact, a necessary condition for the theorem is that the potential 
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levels V's, as well as the kinetic terms R's, be extensive with N, a condition that 
cannot be satisfied by an infinitesimal fraction of particles. 

We believe that our theorem, due to its universality, will have many important 
consequences and we forecast to see in the future many interesting applications. 
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Appendix A. Proof of the Lemma (Eqs. (|53]) - (f5T)) ) 

To prove Eq. let us now come back to the EPR. Recalling that after N t jumps 

the system visits N t + 1 configurations, we can decompose E (.M' '*-*) as 



nr. dCi JO 

- y o <r) ]*» A W... e [<-<]^ A (r) j {A1) 



x el 



where: Q is the set of trajectories extracted with the probability p Nt (the Markov 

?M T/(r) i/to 
l N t ' v o ' ■ ■ ■ i v N t i 



chain starting from n ), and - recalling Defs. (f20|) - Rq P \ . ■ . , R^ , Vq > ■ ■ ■ i ^jv 



(r) (r) 

and , ... ,X N , are the sequences of the the kinetic terms, the potential values, 

and the signs, occurred along the r-th trajectory, respectively, while 6q T \ . . . ,S^'_ t 
are the sequences of the local living times extracted with the density probability 
p(S^ r \ . . . , given by (see note PH) 

p(#>, . . . , 41,) = RPe-<W • • • <_ ie -<-< 5 - x e -<*S2 (A.2) 
and is the residual time given by 



*S=*-W,-,^-i)- ( A -3) 



Finally, as in Sec. Ill, the symbol J" * is a shorthand notation to indicate that the N t 
time intervals Sq T \ . . . , J]^ _ t are constrained by < 5q + . . . + _ x < i. Note that, 

given the intervals Oq , . . . ,ON t -i, the last factor e N * N * in Eq. (|A.2|) represents 
the probability (and not a pdf !) for the system not to jump during the residual time 

interval 5$ (notice the absence of the factor R^) ) . Finally, let rj^ , . . . , 4' be the 
hopping terms occurred along the r-th trajectory. According to the probability of the 
Markov chain, Eq. (fT51 ), the probability p^ is given by 

to J r ) 

* R (r) Rto V ; 

Let now r be some random time at which the system performs some task before 
t, e.g., the first time at which the system enters in some target region of the Fock 
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space. Let us consider the expectation E (jAn'^J- We want now to decompose this 

expectation as we have done in (|A.1|) . However here, unlike Eq. (|A.1[) , we have to 
pay attention to two new points: i) t is random variable taking values in [0, t), ii) by 
definition r is also the time at which a further jump takes place and sends the system 
in some region for the first time before t. Given a trajectory starting from tiq, let us 
indicate with N T the number of jumps performed by the system during the interval 
[0, t). Note that now, the set of the N T time intervals (So, ... , Sn t -i, o~n t ), is a set of 
independent random variables being r = do + ■ ■ ■ + <5jv t a random variable whose pdf 
p(r = x) is given by 

P (t = x) = R%le- R WAlW, (A.5) 

where, given x, 5^ (x) is the residual time as given in Eq. (|A.3j) : 

*s = / *-(^ ) +-+«a_ 1 ). (a.6) 

Note the difference with the previous case in which, given a trajectory with N t jumps, 
the sum of the delta's in Eq. (|A.1[) had to give the constant t. Here we are in 
fact free to choose between one of the two sets of independent random variables to 
be integrated over the interval [0,t): (So, ■ ■ ■ , Sjf T -i, £?0 or (<5o> ■ ■ ■ , $n t -i, t). Wc 
choose the latter, the Jacobian being 1. Furthermore, taking into account that, at 
the jump time x, given R(n x -) (= Rn x ) and i?out(^a;-) (= Rout-,N m ), the system can 
go to any of the R ou t(n x - ) configurations of the reservoir with a total probability 
Rout(n x -) / R(n x -), and by using Eqs. (|A.5[) and (|A.6[) . we arrive at 
r* 



{m^) = / dx £ / c«w . . . (s£\ . . . , 



(r) 



,(r) 



xe L-o--^ J 4 o A ir;... eL «^-^ J0 ^ A ir; ^ t (n x -), (A.7) 

where fi^ is the set of trajectories which, at the instant x, and not before, reach the 
target region. In turn, from Eqs. (|A.2[) and (|A.4[) . we see that the above expression 
can be rewritten as 



5 W WO 



where pi n (^cT^ • ■ ■ ' ^jvi-ij an d PhrN^' are the P^f anc ^ the trajectory's probability, 
in which appear only the kinetic terms (i?i n; o • • ■ , Rin;N w ) and (i?i n ;0 • • • > Rm-N m -i), 

(r) 

respectively. By observing that p in . N is normalized to 1 over the set U x , and by 
comparing with Eq. (jA.lj) we see therefore that Eq. (|A.8j) can be read as 

E (m^) = J a dxE (M%£ a R out (n x -)\ r = x). (A.9) 

The generalization to the case of k random times is trivial. 

Finally, to prove Eq. (|55|) . we observe that: i) the conditional expectation in the 
numerator of the lhs of Eq. |5^|) is - by definition of conditional expectation - built 
with a probability normalized over the space of trajectories leaving the cavity at the 
time x\, ii) as x\ goes to oo, the set of the weights M.^°' Xl ^ (weighted trajectories) not 
included in the expectation of the numerator with respect to those of the denominator, 
approaches 0. The proof of Eq. |F7|) is equal. 



Quantum Phase Transitions induced by Infinite Dilution in the Fock Space 



21 



References 

[1] See for example: S. Sachdev Quantum Phase Transitions (2000) Cambridge University Press. 
[2] Ostilli M and Presilla C, 2006 J. Stat. Mech. P11012 

[3] Beccaria M, Presilla C, De Angclis G F and Jona-Lasinio G, 1999 Europhys. Lett. 48, 243 

[4] Ostilli M and Presilla C, 2005 J. Stat. Mech. P04007 

[5] A. J. Bray and M. A. Moore, J. Phys. C 13, L469 (1980). 

[6] M. Mezard, G. Parisi, and M. A. Virasoro, Spin Glass Theory and Beyond (World Scientific, 
Singapore, 1984). 

[7] T. Obuchi, H. Nishimori, and D. Sherrington, J. Phys. Soc. Jpn. 76 054002 (2007). 
[8] See for example: Shiryayev A N, 1984 Probability (New York: Springer) 

[9] B. Derrida, Phys. Rev. Lett. 45, 79 (1980); D. J. Gross and M. Mezard, Nucl. Phys. B 240, 431 
(1984). 

[10] T. Jorg, F. Krzakala, J. Kurchan, and A. C. Maggs, Phys. Rev. Lett. 101, 147204 (2008). 
[11] Ostilli M and Presilla C, in preparation. 

[12] Godsil, C, Royle, G. Algebraic graph theory, in Graduate Texts in Mathematics, vol. 207 

(Springer- Ver lag, New York 2001). 
[13] See for example: N. Trivedi and D. M. Cepcrley, Phys. Rev. B 41, 4552 (1990). 
[14] Ostilli M and Presilla C, 2004 J. Phys. A 38, 405. 

[15] However, in general, it is not necessary to require that H be the Fock representation of some 
true Hamiltonian operator H. 

[16] We point out that in the proof of the theorem there appears only the invariant measure if 
associated to the cavity space, and not the "flat" measure p. It is possible therefore to 
formulate the theorem without requiring that be pN — > 0, but only ttN — > 0. This more 
abstract formulation is a bit more general and allows to consider systems for which not all the 
kinetic terms K(n) be extensive, but only a finite fraction of them, or even systems having, 
in the Fock space, a finite fraction of kinetic terms going to in the thermodynamic limit. 

[17] Note in general that, the pdf defined as p(s) = exp(— sR)R, is normalized to 1 only in 
the interval [0,+oo), whereas the correct normalized pdf for the finite interval [0,t) is 
p(s) = exp(—sR)R/(l — cxp(— tR)). However, since we are interested only in the limit 
t — > oo, from now on, throughout the paper, to avoid too cumbersome equations, the effective 
approximation R. ~ R./(l — exp(— tR")) will be left as understood. 

[18] Sometimes, to stress the fact that a configuration n belongs to the reservoir or cavity spaces, 
we will make use of the symbol n or n, respectively. 

[19] Notice the null role of the random variables t\, . . . , t^ in the limit t — > oo. 

[20] The precise quantity giving exactly the average number of jumps would imply the knowledge of 
where is the GS energy of the system with V = 0. 

[21] It is possible to reach more formally this result through the decomposition Qt(k) = 
Y2 < N=kQ(k\N)Pt(N), where Pt(N) is the probability to have a total of N jumps during 
the interval [0, t) and Q(k\N) = 'Ppf(k) is the probability for the Markov chain to visits k 
times the cavity. Then, by using for P t {N) ~ ({R)t) N /N\ exp(-(R)t) (see gD]), we see that, 
for t > kt c Q t (k) -> (vr(n)) fc . 

[22] We point out that, in the remind of the proof of the theorem, just to avoid an excessive and 
cumbersome formalism, we will leave the symbols t\ , . . . , tj, as random variables which, for 
the final result, do not play any role. The skeptical reader may rewrite all the equations with 
fixed values ti = j/i, . . . , tj. = y^, and, only at the end of the calculation, to integrate the 
corresponding conditional expectation multiplied by the pdf p(t± = yi, . . . , t& = j/j.). Since the 
conditional expectation, at the leading order, does not depend on the variables j/i, . . . , y^, it is 
clear that, at the leading order, carrying on the calculation with the expectation conditioned 
by the t's, is completely useless. We stress again that, in our decomposition, the role of the 
random variables t's is very different from the role of the random variables r's, and we need 
to use the Lemma only for the latter. 

[23] In fact, the property i) comes immediately from the EPR and the definition of the functional 
A^' '*' which, in the case of the system with Hamiltonian H** , becomes simply A^f**'' '^ = 1. 
Concerning the property ii), since the cavity is connected to the reservoir, there are 
configurations of the cavity such that R* ut > 0, so that, since V* > R" , it must be also 
E* > E" = 0. 

[24] C. Presilla and G. Jona-Lasinio, at the time in which the EPR was derived, emphasized 
the exactness of the EPR, and hypothesized possible advantages not featured by other 
approximate representations. 



